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We present thermal and electrical transport measurements of low-density (10 14 m -2 ), mesoscopic 
two-dimensional electron systems (2DESs) in GaAs/AlGaAs heterostructures at sub-Kelvin temper- 
atures. We find that even in the supposedly strongly localized regime, where the electrical resistivity 
of the system is two orders of magnitude greater than the quantum of resistance /i/e 2 , the ther- 
mopower decreases linearly with temperature indicating metallicity Remarkably, the magnitude of 
the thermopower exceeds the predicted value in noninteracting metallic 2DESs at similar carrier 
densities by over two orders of magnitude. Our results indicate a new quantum state and possibly 
a novel class of itinerant quasiparticles in dilute 2DESs at low temperatures where the Coulomb 
interaction plays a pivotal role. 



I. INTRODUCTION 

The thermopower or Seebeck coefficient S of a system 
is the electric voltage Vth generated in response to an im- 
posed temperature difference AT across its ends. In an 
electron system thermally connected to a phonon bath, 
the primary contributions to V t h are from "phonon drag," 
which is driven by electron-phonon scattering, and the 
diffusive kinetics of the electrons themselves in order to 
maintain local thermal equilibrium. The latter is domi- 
nant at sufficiently low temperatures when most phonons 
freeze out and is called the diffusion thermopower {So)- 
Within the semiclassical Boltzmann framework, the dif- 
fusion thermopower is given by the Mott relation pQ and 
connects it to the conductivity a of the system as 



S d = 



Vth 
AT 



3q 



d\no 
dE 



(1) 



where ks is the Boltzmann constant, T is the average 
electron temperature, q is the charge of the carriers, E 
is the energy, and /i is the chemical potential of the sys- 
tem. Thus, S is sensitive to the energy- derivative of a 
and consequently various system parameters such as the 
electronic density of states (DOS) and the momentum 
relaxation time r. This sets it apart from the resistivity 
as a spectroscopic tool and the physics it probes. For 
instance, the T dependence of S is markedly different for 
metals, Anderson insulators, and gapped insulators: For 
metals, where there is a continuous DOS and free charge 
carriers, S(T) decreases to zero linearly as T — >• 0; for 2D 
Anderson insulators which are characterized by a sharp 
mobility edge separating the conducting (extended) and 
nonconducting (localized) states in an otherwise contin- 
uous DOS, S(T) varies as T 1 / 3 ; in the Efros-Shklovskii 
regime, where Coulomb interactions cause the opening 



of a soft gap at the Fermi energy in the DOS, S(T) — >• 
constant as T — > 0; and finally, for hard gapped insula- 
tors S(T) diverges as 1/T. These distinctions are far less 
pronounced in the resistivity p(T). Another virtue of S 
is its equivalence to the entropy per carrier, thus provid- 
ing a probe of this abstract quantity. We note here that 
the Mott formula [Eq. ([I])] was originally derived for a 
strongly degenerate, noninteracting metal and hence is 
not expected to be valid when strong interactions are 
present. 

The theromopower has proved to be a powerful tool to 
probe many-body phases in strongly interacting 2DESs, 
particularly those involving transitions from insulating to 
either metallic or fractional quantum Hall liquid phases. 
In the presence of strong transverse magnetic fields, the 
thermopower has been employed to investigate the nature 
of the insulating state near the v = 3/2 fractional quan- 
tum Hall state [2], the existence of an energy gap at the 
v = 5/2 fractional quantum Hall state [3], and the pos- 
sibility of a Wigner crystalline ground state near v = 1 
in a bilayer GaAs/AlGaAs hole system [4]. The possi- 
bility of a metallic ground state in low-density 2DESs 
at zero magnetic field, where the metallicity is driven 
purely by Coulomb interactions, has however been far 
more controversial [5 j. Several attempts to identify such 
a metal have been made with thermopower as a probe in 
the apparently metallic regime, notably in low-density 2D 
electron or hole systems [6j [7] and high- mobility Si MOS- 
FETs [8]. While these studies indicate a definite change 
in the conduction mechanism and/or critical behavior in 
thermopower near the transition point between insulator 
and apparent metal, an unambiguous demonstration of 
a metallic state has never been achieved. Much of the 
uncertainty could be due to the very nature of the tran- 
sition which has been suggested to be an inhomogeneity- 
driven, classical percolation transition [9j [10] rather than 
an interaction-induced one. 
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heterostructures in which the distance between the 
dopant layer and 2DES is 40 nm. The as-grown mobil- 
ity was 220 m 2 /Vs at carrier densities of 2.2 x 10 15 m~ 2 . 
The mesoscopic 2DES is defined using a top-gate of di- 
mensions L x W = 2/im x 8/im, which enables us to 
tune its density n s by applying a gate voltage V g . The 
n s -V g calibration is obtained via an edge state reflection 
technique [17] • We measure the 2DES resistivity p2DES 
in a 4-probe setup by passing a small excitation current 
I ex = 100 pA at / = 7 Hz and detecting the output V ex 
using a lock-in amplifier. To measure S we impose a tem- 
perature difference AT across the device by means of a 
heating current 1^ = 4-5 fiA at fh = 11 Hz and detect 
the thermovoltage Vth at 2/^. AT across the device is 
measured using large (~ 10/i m) b ar gates following ref- 
erences [18] and [19] (see Sec. II A for further details). In 
the results reported here AT never exceeds 20 mK (see 
Fig. [3]). In our experiments P2DES-, Vth and AT are all 
measured simultaneously. We perform measurements on 
three devices and observe similar results in all. 



FIG. 1: (Color online) (a) False-color SEM image of a typ- 
ical device, (b) Schematic of the device layout. D denotes 
the mesoscopic device, BG1-4 are macroscopic bar gates used 



for thermometry (see Sec. II A), and 1-8 represent ohmic con- 
tacts. A heating current Ih establishes a temperature gradient 
along the length of the device. 



Recently, the possibility of a metallic phase in strongly 
interacting 2DESs in zero magnetic field has resurfaced 
through transport experiments with micron-scale 2DESs 
in GaAs/AlGaAs heterostructures [TTUT4] . The motiva- 
tion behind using micron-sized or "mesoscopic" samples 
was to circumvent the influence of the long-ranged dis- 
order known to exist in molecular-beam-epitaxy-grown 
GaAs/AlGaAs heterostructures [15J[T6]. This greatly re- 
duces the tendency of the 2DES to fragment into puddles 
and, consequently, the chances that the transport be gov- 
erned by percolation through these puddles. Indeed, at 
low n S: Baenninger et al [13] observed striking behavior 
in mesoscopic 2DESs: The T dependence of the resis- 
tance weakened drastically either saturating to a finite 
value or decreasing with decreasing T even though the 
absolute value of resistivity exceeded h/e 2 by several or- 
ders of magnitude. The aim of the present work is to 
investigate the nature of the DOS at the Fermi energy 
and many-body phenomena in these systems using ther- 
mopower measurements. 
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II. EXPERIMENT 



Figure [TJa) shows a micrograph of a typical device. 
The devices are fabricated in Si-doped GaAs/AlGaAs 



FIG. 2: (Color online) Density and temperature dependence, 
respectively, of S in the high-n s regime. The dashed lines 
show S d [Eq. S] 
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A. Thermometry using bar gates 

We measure the local electron temperature by measur- 
ing the thermovoltage between large 2DESs 10 /am x 
30 jam) defined by metallic gates. The layout of the bar 
gates can be seen in Figs.[]Ja) and[]Jb). At high n s when 
interaction or localization effects are negligible and the 
2DES is well described as a noninteracting, Drude-like 
metal, Eq. ([!]) reduces to, 



s d 



■nk 2 B Tm 1 
3eh 2 



a 



(2) 



Here ks is Boltzmann's constant, T is the average tem- 
perature = (T e +Tl)/2 with T e and Tl the electron and 
lattice temperatures, respectively, m is the effective elec- 
tron mass in GaAs , — e is the electronic charge, h is 
Planck's constant divided by 27r, n is the 2DEG den- 
sity, and a = where r is the momentum relax- 
ation time. To measure ATi, say [see Fig. [TJd], where 
the local electron temperature is T e i, we differentially 
bias BGi and BG2 and detect the the signal Vati be- 
tween contacts 2 and 7 at twice the heating frequency. 
The difference in thermopowers between BGi and BG2, 
AS d = V^/AT = V ATl /(T el — Tl). Substituting the 
expression for S d from Eq. ^ we obtain 



^el = ^ATx 



3eH 2 



7rkgm(l + a) \ni ri2 



(3) 



where n\ and ri2 are the 2DEG densities beneath BGi 
and BG2, respectively, and are obtained via an edge-state 
reflection technique [IT] . Tl is obtained from a RU2O 
thermometer attached close to the sample. A similar 
procedure yields the value of T e 2 from which AT, the 
temperature difference across D, the mesoscopic device 
under study, is obtained as T e2 — T e \ . 

In order to calibrate our setup we have first carried out 
detailed thermopower measurements when D is tuned to 
relatively high n s where p2DES < h/e 2 and interaction 
or localization effects are negligible. Here the 2DES is 
well described as a noninteracting, Drude-like metal and 
its thermopower given by Eq. ([2|. Figure [2^a) shows the 
dependence of S on n s in the regime n s > 3 x 10 14 m -2 , 
measured at base temperature, T = 0.28 K. The dashed 
line represents Sd for a = 1. In Fig. |5|b), we present 
the T dependence of S at two values of n s [= 4 and 
5 x 10 14 m -2 , denoted by the vertical dashed lines in 
Fig. pja)]. At both n s we find S to be well described by 
Eq. pi) . Apart from a consistency check for the measure- 
ment process, Figs. [2|a) and [2^b) also confirm that S is 
dominated by the diffusive component over this range of 
T (< 1.5 K), and that the phonon drag contribution to 
S is negligible. This is not surprising as the hot-electron 
technique greatly reduces the phonon drag contribution 
as was demonstrated by Chickering et al. in Ref. [19]. 



A second validation of our measurement technique is 
shown in Fig. [3^ a) . We see that AT shows two distinct 
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FIG. 3: (Color online) (a) The figure shows the variation of 
AT with the lattice temperature Tl . Also shown is the qual- 
itative dependence of AT expected due to heat dissipation 
by piezoelectric electron-phonon scattering and to the ohmic 
contacts, (b) The figure shows that as In is increased AT in- 
creases approximately quadratically before nonlinearities ap- 
pear at large Ih- 



dependencies on T: A slower decay when T < 0.6K and 
a faster decay when T > 0.6K. The behavior of AT 
can be accounted for simply by a combination of power 
dissipated to the ohmic contacts, Po? an d through piezo- 
electric electron-phonon coupling, P ep . The net power 
dissipated P = Po + P ep = I^R- Here Ih is the (con- 
stant) heating current and R the resistance of the heating 
channel [shown as red arrow in Fig.^b)]. In the range 
of temperatures explored, R is completely dominated by 
impurity scattering and hence shows no temperature de- 
pendence. Consequently, P is independent of T. Assum- 
ing P = K^Tl-T 2 ) and P ep = K 2 (T e 5 -T|) [20^ where 
Ki and K 2 are system-dependent constants, we obtain 
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AT = T e - T L <x (^(Te + T L ) + K 2 {Tt + T 3 e T L 



T 2 T 2 



T e Tl 



n 



))- 



(4) 



The qualitative behavior described by Eq. Q is shown 
as a broken line in Fig. [3^a) and is seen to adequately 
describe the observed dependence. We mention here that 
even in the absence of a voltage on the gate finger D in 
Fig.[TJb), ATi = 0; i.e., the electrons relax to Tl over the 
distance between the heating channel and cold end of the 
device. Consequently, AT = T e 2 — Tl does not change 
as the 2DES is tuned between the high- and low-density 
regimes. A third check is shown in Fig. |3jb) where we 
plot AT = T c — Tl against the heating current 1^. We see 
that as expected, AT grows as « l\ before saturating, 
presumably due to nonlinearities arising from the large 
values of Ih- 



III. RESULTS 



drag to set in at lower T. 

Figure [5|a) shows the n s dependence of S at T = 
0.28 K. There are two salient features which need to 
be noted: First, while S increases with decreasing n s , 
the decrease is much stronger than the ~ l/n s behav- 
ior expected from Eq. (J2|; S increases to about two or- 
ders of magnitude above the Mott value at the lowest 
n s . Second, the increase in S is oscillatory rather than 
smooth, even changing sign occasionally, although the de- 
tailed nature of oscillations is highly device dependent. 
Coulomb blockade due to electron puddles in an inho- 
mogeneous charge distribution could lead to oscillatory 
thermopower [25 , but such a scenario is unlikely because 
(1) in spite of occasional change in sign, the oscillations 
are primarily one sided (negative) and very asymmetric 



As n s in the mesoscopic region is reduced by mak- 
ing the top-gate voltage more negative, p2DES increases 
rapidly below n s ~ 2 - 3 x 10 14 m~ 2 (inset Fig. [5] and 
Refs. [13 and [14]). The T and n s dependencies of 
P2DES m such low-n s , mesoscopic 2DESs are very differ- 
ent from those of their macroscopic counterparts [I3j [14] . 
In Fig. 4^a) we see, similar to what was observed in 
Ref. [13 . that at low n s values, where the 2DES is ex- 
pected to be localized, p2DEs(T) does not diverge as T — > 
but saturates at values ^> h/e 2 below « 0.8 K, indi- 
cating noninsulating behavior. The precise mechanism 
for resistivity saturation is debated, ranging from quan- 
tum tunneling between multiple electron puddles [2TJ [22] 
to defect migration in a spontaneously broken symmetry 
phase [13] [14] . Figure pFb) shows S as a function of T at 
the same n s as Fig.[4^af Its behavior is, however, generic 
to all n s in the "localized" regime: At n s < 2 x 10 14 m -2 , 
we find that S always decreases with decreasing T despite 
the fact that p2DES ^ h/e 2 . The nature of its decrease 
depends on the range: At T < 0.8 K, S varies linearly 
with T, as in the case of a metal. This linear dependence 
is observed down to the lowest measured n s where the 
electrical resistivity can be as high as 200 — 300 x h/e 2 . 
When T > 0.8 K, S is seen to grow at a much faster 
rate, S ~ T 4 , as a result of which S reaches extremely 
large values ~100 mV/K at 1.4 K. The T 4 dependence 
of S in GaAs/AlGaAs hetero junctions at low T is usu- 
ally understood in terms of phonon drag [4j [23] , but its 
absence at the same T at high n s in our devices (see 
Fig. [2| suggests that the present instance warrants closer 
inspection. We note, however, that screening by free car- 
riers in GaAs/AlGaAs heterostructures can significantly 
alter the polaron (electron-LO phonon) binding energy 
and scattering rate [24], and it is possible that the ab- 
sence of free electrons in the low-n s regime allows phonon 



n s (10 1 V 2 ) 




FIG. 4: (Color online) (a) In the highly resistive regime P2DES 
becomes almost T independent below ~ 0.8 K. (b) S in- 
creases linearly below T < 0.8 K beyond which it increases 
as T 4 , in agreement with earlier reports of thermopower in 
GaAs/AlGaAs heterostructures (Ref. [4]). 
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n s (10 15 m" 2 ) T(K) 



FIG. 5: (Color online) (a) S vs n s for 0.28 K < T < 0.7 K. The broken green line shows S d [Eq. ^} at 0.28 K. Inset: p 2 DES vs 
n s at the same T values; there is little T dependence in this range, (b) Low-T linear variation of S. Inset: Descriptions based 
on variable-ranged hopping, where S is expected to decay to zero as T 1 / 3 , do not adequately describe the observed data. 



around zero, and (2) none of this structure is seen in 
P2DES [shown in inset to Fig. |5ja)] which grows mono- 
tonically and is, by and large, featureless. In fact, the 
energy derivative of p2DEs( n s) fails to account for the 
oscillations (see Sec. |III A| both qualitatively and in mag- 
nitude, which indicates a breakdown of the semiclassical 
Mott picture [Eq. Q]. While the precise origin of the os- 
cillations remains unexplained, a many-body interaction 
effect, for example a reentrant order-disorder transition 
in a broken-symmetry phase, cannot be ruled out [26 . 

Figure [5jb) embodies the key result of this work. We 
inspect the linear regime of S(T) and show linear fits to 
the data that pass through the origin. This suggests the 
absence of any gap, hard or Efros-Shklovskii type, at the 
Fermi energy. In the inset we see that the same data 
plotted against T 1 / 3 cannot be described as a line going 
through the origin, showing clearly that the observed be- 
havior is not consistent with Mott-type variable-ranged 
hopping in an Anderson insulator. We emphasize the 
difference between our result and earlier thermopower 
results near the apparent metal-insulator transition in 
macroscopic 2D systems: Apart from this being the first 
clear observation of metal-like thermopower in the low- 
n s regime, we find this behavior to exist over an ex- 
tended range of n s and p2DES (up to ~ 300 x h/e 2 
which is the upper limit of our measurements), unlike 
previous studies which explore only the critical regime 
(p2DES ~ 0.1 - 1 x h/e 2 ). 



related by the Mott formula [Eq. ([TJ]. To make a quan- 
titative comparison with the Mott result applied to our 
system, we evaluate the energy- derivative of In a using 
the Hartree expression for the energy of the 2DES 



E = Ek + Eq 



ma b 



(5) 



where as is the effective Bohr radius in GaAs ^11 nm. 
Using Eq. ^ the expression for Smott reduces to: 



WT = ^(l+r s /2)- 1<nnp 



3eh 2 
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A. Breakdown of the Mott relation 

In this section we show that the electrical and thermal 
transport in the low-n s , mesoscopic 2DES are not simply 



FIG. 6: (Color online) The figure shows S and 
Smott [Eq. |l])] at T = 0.28 K. The S trace is also shown in 
Fig. [5ja). Note that the left and right vertical axes differ by 
a factor of 10. 
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n s (10 15 m" 2 ) 

FIG. 7: (Color online) The figure shows S and p2DES over 
the same n s range when p2DES ~ h/e 2 . The oscillations in S 
are strikingly present even at these low resistivities. 



In Fig. [6] we compare the measured thermopower S 
at T = 0.28 K to that expected from the Mott formula 
Smott [the former is also presented in Fig. |5ja)]. We 
note that S and Smott disagree qualitatively and quan- 
titatively. First, S is nearly two orders of magnitude 
larger than Smott and second, though Smott captures 
some of the broad features in 5, there is additional struc- 
ture in the latter. The same is true for both of the other 
devices measured. 

In Fig. (7) we see that the oscillations in S are observed 
even when p2DES ~ h/e 2 . The strong oscillations and 
even sign changes in S are unaccompanied by any corre- 
sponding structure in p2DES- This strongly suggests that 
these oscillations are not a result of Coulomb blockade in 
electron puddles. 



IV. DISCUSSION 

Two outstanding questions remain: First, can the na- 
ture of the metallic state be understood within an effec- 
tive semiclassical Boltzmann framework? And second, 
why is S two hundred times larger than that expected 



from the non-interacting model at low n s ? In the Drude- 
like metal described by Eq. (|2|, the resistivity and ther- 
mopower are both inversely proportional to the density 
of delocalized quasiparticles which carry heat and elec- 
tricity, i.e., S ~ l/n ex ~ P2DES, where n ex is density 
of delocalized quasiparticle excitations. The deviation of 
the observed S from the Mott expectation clearly im- 
plies that n ex is not given simply by n s . However, to see 
whether a Drude-like description is valid we have plot- 
ted S as a function of the corresponding p2DES for three 
different devices at base temperature in Fig. [8j Remark- 
ably, in spite of the superposed oscillatory structure, we 
find S to be linearly proportional to p2DES over nearly 
three decades in S in all three devices, providing strong 
evidence of a Drude metal-like character, but with a dif- 
ferent nature of itinerant quasiparticles. The number of 
such quasiparticles can be significantly smaller than n s 
which can readily explain the large magnitude of both S 
and p2DES- Note that Fig. [8] provides further evidence 
against percolative transport via tunneling between elec- 
tron puddles that act as quantum dots since, in such 
cases, S is expected to be proportional to the energy 
derivative of p2DES rather than p2DES itself [25] . 

In the presence of strong Coulomb interactions, en- 
hanced screening and the formation of many-body states 
have been suggested to lead to extended, or at least rela- 
tively less localized, electron wave functions in the ground 
state thereby modifying the single-electron conduction 
mechanism (see for example Ref. [27]). Moreover, a com- 
mon effect of interactions at low temperatures is to in- 
duce broken-symmetry states such as Wigner crystals, 
stripe/bubble phases [28 etc. Recently, the metal-like 
electrical transport observed in Refs. 13 and Koushik 
was attributed to the melting of spontaneous symmetry- 
broken states [14^ , where the melting transition proceeds 
through the proliferation of topological defects that form 
delocalized quasiparticles at sufficiently low temperatures 
where zero-point fluctuations dominate [29]. The signa- 
ture of the Berezinskii-Kosterlitz-Thouless (BKT) melt- 
ing transition was indeed observed in all our devices [in- 
sets of Figs. [8](a) -|8jc)], where the resistivity p2DES 
varies as p2DES °c exp(— A/^n s — n c ), A being a con- 
stant and n c the melting density [14 . The correspond- 
ing delocalized excitations can be responsible for thermal 
transport as well, and being much smaller than the num- 
ber of electrons in the system, may lead to very large 
thermopower. We note that though Wigner crystallisa- 
tion is expected at much lower n s in disorder-free sys- 
tems [30], the presence of disorder may "pin" the order 
and allow finite grain-size Wigner crystallites at much 
higher n s [3T] . 

Thus, in conclusion, the effects of the many-body 
Coulomb potential appear to manifest in a novel metal- 
lic phase in dilute mesoscopic 2DESs, with giant ther- 
mopower that decreases linearly as T — ^ 0. While we 
cannot rule out the possibility that the system turns in- 
sulating at even lower T, the fact that it persists down to 
lowest experimental n s explored here indicates that the 
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FIG. 8: (Color online) (a)-(c) In all three devices measured there is a clear linear envelope between S(n s ) and P2des(ji s ) as is 
expected in simple Drude-like metallic system. However, the n s dependence of neither obeys a simple Drude description. The 
insets show that p2DES in each instance is well described by a Berezinskii-Kosterlitz-Thouless model (gray lines) suggesting 
that topological defects of an underlying ordered 2DES could be the mediators of transport in the dilute, mesoscopic 2DES. 



metallic phase might be the true ground state of strongly 
interacting 2DESs. 
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